Wigner Glass, Spin-liquids, and the Metal-Insulator Transition 
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Recent experiments on the two dimensional 
{d — 2) electron gas in various semiconduc- 
tor devices have revealed an unexpected metal- 
insulator transition and have challenged the 
previously held assumption that there is no such 
transition in two dimensions. While the exper- 
iments are still at the stage of rapid develop- 
ment, it is becoming evident that they cannot be 
understood from the conventional perspective of 
weak interactions. In the present paper, we pro- 
pose the following. (1) The low-density insulat- 
ing state is the Wigner Glass, a phase with quasi- 
long-range translational order and competing fer- 
romagnetic and antiferromagnetic spin-exchange 
interactions. (2) The transition is the melting 
of this Wigner Glass, disorder being the agent 
allowing the transition to be second order. (3) 
Within the Wigner Glass phase, there are at least 
two, distinct magnetic ground-states, a ferromag- 
netic state at very low electron density and a spin- 
liquid state with a spin pseudo-gap at higher den- 
sities. (4) The metallic side of the transition is a 
non-Fermi liquid. These conclusions are encap- 
sulated in Figure 1 which presents the proposed 
phase diagram as a function of disorder strength 
and density; we also suggest experimental signa- 
tures of the various phases and transitions. 

Although there is not complete agreement between ex- 
periments, a number of striking features have emerged. 
Most importantly, the temperature dependence of the re- 
sistivity changes at a critical carrier concentration (ric): 
for lower carrier concentration {n < ric), it is insulat- 
ing, implying a diverging resistivity as the temperature, 
T ^ 0; for higher carrier concentration (n > ric), it is 
metallic, the resistivity decreasing with decreasing tem- 
perature. Therefore, ric appears to signify a zero temper- 
ature quantum phase transition between two fundamen- 
tally distinct states of matter. It was previously thought 
that such a metal-insulator transition is not possible in 
two dimensions |^]. Indeed, previous experiments, car- 
ried out over the past two decades on systems with higher 
electron density where the effects of electron-electron in- 
teractions are significantly weaker, strongly supported 
the idea that a metallic state in two dimensions is not 
possible. 

In some of the recent experiments, a striking reflec- 



tion symmetry is observed It is striking because, 
a priori, there is no particular symmetry between the 
two distinct states of matter, the metal and the insula- 
tor. Over a significant range of n around Uc, one finds 
that g*{Sn,T) = l/g*{-Sn,T), where g* is the conduc- 
tance scaled by the value at the critical concentration, 
which is of order e^/h, e being the electronic charge and 
h the Planck's constant, and (5„ — {n ~ Uc) / Uc- Concomi- 
tantly, the data obeys scaling, that is, the data can be 
collapsed onto two universal curves, one for the metallic 
samples and one for the insulating samples, when it is 
plotted against the scaling variable [To((5„)/T]^/^, where 
X is a critical exponent, and the scaling parameter is 
To{Sn) ^ \Sn\^- The resistivity as a function of the elec- 
tric field shows a similar scaling behavior Q, and the 
current-voltage characteristic shows equally striking re- 
flection symmetry; the data on the insulating side can be 
mapped on to the metallic side by a simple reflection. 

The immediate conclusion one can draw from these ex- 
periments is that the scaling and reflection properties are 
due to a correlation length diverging at a quantum crit- 
ical point ^ corresponding to a quantum phase transi- 
tion. Moreover, the observation that the range of carrier 
concentration over which reflection symmetry is obeyed 
shrinks with decreasing temperature is consistent with 
a zero temperature quantum critical point There- 
fore, although the phenomenological scaling theory ||] 
appears to be valid, the following key questions remain 
unaddressed: What is the origin of the quantum critical 
point! What is the nature of insulating and conducting 
phases'? In this paper we substantively address the ori- 
gin of the quantum critical point and the nature of the 
insulating phases, and offer a conjecture regarding the 
conducting phase. 

There have been theoretical attempts to explain the 
experiments in terms of a conventional Fermi liquid ap- 
proach ||l^] pioneered by Finkel'stein At best, such 
a perturbative renomalization group approach can pro- 
vide hints concerning the existence of a quantum critical 
point, since this point lies outside the perturbative regime 
]l2t . In addition, the nature of the strong coupling phases 
discussed below are such that they simply cannot be ac- 
cessed in perturbation theory. Because the experiments 
do not indicate a finite superconducting transition tem- 
perature in samples with n > ric, it remains to be seen 
if the "metallic phase" is actually a superconductor in 
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disguise (T3-15|. The sensitivity of the resistivity data to 
magnetic fields at the temperatures probed in these ex- 
periments [pljl^] appears inconsistent with any classical 
explanation. 

New insight into this problem of the disordered two di- 
mensional electron gas can be obtained if it is approached 
from the strong interaction (low density) limit the 
requisite quantum critical point is due to the zero tem- 
perature continuous melting of a new state of matter, 
called the Wigner glass, which, as shown below, is dis- 
tinct from either an Anderson insulator, or a Mott in- 
sulator |5j. At this quantum critical point scaling must 
hold. In contrast, in the absence of disorder, the in- 
sulating triangular Wigner crystal, which minimizes the 
Coulomb energy, quantum melts via a first order transi- 
tion at which scaling does not apply. Thus, the situation 
is radically different with and without disorder. It is the 
principal purpose of this paper to analyze the nature and 
the consequences of the existence of the Wigner glass. 

The presence of electronic spin degrees of freedom dis- 
tinguishes the Wigner glass from the similar Bragg glass 
psf which has been discussed in other contexts. The 
spin degrees of freedom lead to a magnetic Hamiltonian 
which is highly frustrated due to the presence of signif- 
icant ring-exchange processes. Frustration, in turn, can 
lead to novel magnetic phases, including a "spin-liquid" 
| p9[ , where a "spin-liquid" is a state in which quantum 
fluctuations are so strong that even at T = no sym- 
metries are broken. As the Wigner glass melts, the elec- 
trons are no longer power-law ordered at the lattice sites, 
but the local crystalline character and the short distance 
antiferromagnetic spin-singlet correlations should survive 
within the scale of the correlation length, which is large 
close to the quantum critical point. Most likely such a 
state cannot be a Fermi liquid due to its strong singlet 
correlations. Moreover, the observed conducting state 
cannot be a Fermi liquid for it would localize if it were 
A simple model of a non- Fermi liquid has a per- 
fect metal ground state if the interactions arc sufhciently 
strong po|] , and the interactions are certainly strong in 
the present experiments. We conjecture that the con- 
ducting state is a non-Fermi liquid metal with strong 
spin singlet correlations. 

For a Wigner glass, in the presence of disorder, it fol- 
lows from the work of Giamarchi and Le Doussal fl^ ] 
that the equal-time correlation function 



Cg{R) = (e*G-"(R-)e-'G-u(o)) _ _^ R^oo, (1) 

where u(R, t) is the displacement at the lattice site R, 
and 77 is a critical exponent. The angular brackets denote 
the groundstate average and the overline the disorder av- 
erage. Therefore, the system exhibits quasi-long-range 
crystalline order characterized by power law Bragg peaks 
at the reciprocal lattice vectors, G, of the triangular lat- 
tice, instead of the (5-function peaks of a crystal with true 



long-range order. 

This quasi-long-range crystalline order is possible only 
if unbound dislocations are not generated by disor- 
der. Giamarchi and Le Doussal have given a self- 
consistent argument that this is indeed the case for a 
range of parameters, but, in d = 2, this is still an open 
question; for d = 3, this idea has received further sup- 
port, recently [^,^. Pending further investigation, we 
shall assume the existence of the Bragg glass phase in 
d = 2. Even if the self-consistent argument of Ref. p^ ] 
turns out to be incorrect, dislocations are expected to 
be exponentially rare and so of not great practical 
importance. On the other hand, at very low carrier con- 
centrations, corresponding to very high disorder, disloca- 
tions will proliferate, destroying this phase, and the re- 
sulting state is indistinguishable from the Mott- Anderson 
insulating phase. Thus, the Wigner glass should exist 
only at intermediate carrier concentrations. It is also 
worth noting that the dissipative properties of this glass 
at the melting point should be very different from the 
melting of the Bragg glass because of the continuum of 
low lying fermionic excitations. 

A crucial question is whether or not the Wigner glass is 
a distinct state of matter. The findings of Giamarchi and 
Le Doussal imply that it is. Specifically, it is character- 
ized by a power-law crystalline order (See, Eq. (|l])) and 
by broken replica symmetry |^^, as in a spin glass. As a 
consequence, there are many low lying metastable states, 
with barriers between them diverging as a function of the 
"separation" in the phase space. Such a state cannot be 
accessed by perturbative renormalization group methods 

0- . . . 

As the density is increased, the Wigner glass will quan- 
tum melt in a continuous manner, in contrast to the first 
order melting of a Wigner crystal; in the presence of dis- 
order, coexistence of phases is not possible, and, there- 
fore, a first order transition is forbidden ||24| . However, 
because replica symmetry is restored in the melted liq- 
uid, there must still be a phase transition and hence a 
quantum critical point! 

The pure Wigner crystal transition in d = 2 is believed 
to occur at w 37 p5|] , where rg is a measure of the im- 
portance of electron-electron interactions and is defined 
as the average interparticle separation in units of the ef- 
fective Bohr radius as — h^e/mbe^; e is the dielectric 
constant and mi, is the effective mass of the charge carri- 
ers. The situation changes dramatically in the presence 
of disorder, because the elastic Wigner glass can take ad- 
vantage of pinning and lower its energy relative to the 
liquid phase. There are numerical results that indicate 
that the critical value of can be lowered substantially, 
for example, to 7.5 from 37 in the presence of disorder 
pGf . The critical is, of course, a nonuniversal quan- 
tity and can depend strongly on the precise microscopic 
model. 

The local crystalline order of a Wigner glass allows us 
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to estimate exchange energies. For this purpose, it is 
legitimate to treat the glass as a triangular crystal, rec- 
ognizing that the exchange constants will be randomly 
distributed. A general theorem shows that exchanges 
involving an even number of Fermions are antiferromag- 
netic, while those involving an odd number of Fermions 
are ferromagnetic p7| , p8[ . There are estimates of these 
multiparticle exchanges in both solid '^He and the Wigner 
crystal ||2^. The version of the many dimensional WKB 
pO| formalism used here to calculate them is different 
in detail from that of Ref. . The actual potential at 
short distances is, of course, softer than the Coulomb po- 
tential, e^/er, due to the spread of the wavefunctions of 
the electrons perpendicular to the plane. However, we 
have explicitly verified that a more realistic softer poten- 
tial makes a negligible difference (a little larger exchange 
energies) as long as the spread in the perpendicular di- 
rection is less than half the lattice spacing, which is the 
case by a safe margin. This implies that our results are 
independent of the short distance microscopic details and 
depend only on the nature of the potential on the scale 
of the lattice spacing, which typically is of order 300A. 

The semi-classical expression for a p-particle exchange 
energy Jp is 

Jp = Ap?i^o(|^)'^^-^''/^ (2) 

Here Bp is the value of the Euclidean action along 
the minimal action path that exchanges p electrons on 
a triangular lattice, and too is the attempt frequency, 
which can be estimated from the phonon spectrum of 
the Wigner lattice |^ . We ascribe the zero point energy 
of the phonons to an effective oscillator, htoo = 
in units of the Rydberg constant, which is e^/2eaB- The 
prefactor A, in one dimensional WKB problems, can be 
proven to be greater than unity | |3^ and is often quite 
a bit larger than unity; this should hold more generally. 
We set Ap — I, which is an underestimate. Equation (||) 
can be trusted as long as Bp/Ti is much larger than unity. 

We have calculated Bp/h for a large number of ex- 
change configurations (upto p — 18) and have also made 
asymptotic estimates for p ^ 1. Here, we report re- 
sults for only the most compact 2, 3, 4, 5, and 6- 
particle exchanges for the Coulomb potential, of which 5 
and 6-particle exchanges were not previously calculated. 
This required minimization over typically 500-1000 vari- 
ables to obtain satisfactory convergence. If we define 
Bp/n = bprl^^, then 62 = 1-66, &3 = 1-52, 64 = 1-67, 
&5 = 1.91, and = 1.77. The values of bp for p = 2,3, 
and 4 are typically in the range 10-15% lower than the 
previously known estimates [ p9[ ; this is because full min- 
imization was not previously carried out. 

If we retain for simplicity only exchanges up to 4 parti- 
cles, it is easy to show in terms of the spin- 1/2 operators. 



S, that for the perfect triangular lattice, the spin Hamil- 
tonian is 

NN 

H = est. -I- JpfN Si • Sj 

(u> 

NNN Rhombus 

+ Jnnn ^ Si ■ Sj + 4 J4 ^ Gijki, (3) 

{ij) {ijkl) 

where djki = {Si ■ Sj){Sk ■ Si) + {Si ■ Si){Sj ■ Sk) - {Si ■ 
Sk){Sj ■ Si). The first sum is over all distinct nearest 
neighbors, where Jnn — 4(J2 + |J4 — J3). The second 
sum is over all distinct next nearest neighbors, where 
Jnnn — J a- The third sum is over all distinct rhom- 
buses. Jnn is antiferromagnetic for < 39, and fer- 
romagnetic for larger Tg. Its magnitude, using e = 7.7 
and mb = 0.2 (appropriate for Si-MOSFET) is 5.2 K, 1.6 
K, 0.4 K, 0.06 K for = 5, 7, 10, and 15, respectively. 
Note, however, that while Jnnn / Jnn ~ 0.4 at rs = 5, 
the ratio approaches 1 at about rs = 29. This is a highly 
frustrated antiferromagnetic spin Hamiltonian. Numer- 
ical evidence |^ suggests that the ground state of this 
Hamiltonian is "quantum disordered" for rg ^ 39, where 
"quantum disordered" means that spin-rotational sym- 
metry is not broken, even at T = 0. Because the actual 
Wigner glass is not translationally invariant, any quan- 
tum disordered state is a spin liquid; various dimerized 
states, which break translational symmetry, are not dis- 
tinct states of matter. However, because of disorder, the 
spin liquid only exhibits a spin pseudo-gap state, rather 
than the spin gap state of the pure system . Other 
magnetic structures may develop for ^ 39, but they 
are likely to be predominantly ferromagnetic. 

In the current experiments cited earlier, the critical 
values of rs at which the metal-insulator transition takes 
place are in the range of 5-15. (The character of this 
transition will be very different if it takes place for larger 
values of r^ for which exchange is of insignificant impor- 
tance at experimentally realistic temperatures.) Because 
the exchange energies at rg ~ 10 are about 0.4 K, the 
spin-singlets can be broken with modest magnetic fields 
of order IT, and the ground state should be very sen- 
sitive to such fields. Once the spin-singlets are broken, 
the elctrons can be described as a Fermi liquid, which 
should localize in d = 2 due to disorder. Indeed, it has 
been observed that the metallic behavior in Si MOSFETs 
can be suppressed by such fields In the absence of 
the magnetic field, the continuous nature of the melting 
of the Wigner glass allows the spin correlations on the 
scale of the correlation length (large close to the melting 
transition) to be similar to those in the Wigner glass. In 
particular, the spin pseudogap character should survive 
the melting process, and the drop of the resistivity as 
a function of temperature seen in experiments is due to 
the reduction of the scattering rate of the carriers as the 
temperature is lowered due to the opening of the spin 
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pseudogap. 

The quahtative global phase diagram at T = is 
sketched in Fig. |^ with D as the measure of the strength 
of the quenched disorder. 



AMORPHOUS INSULATOR 
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FIG. 1. The T = phase diagram as a function of disorder 
D and Ts- The point X is a tetracritical point. The arrowed 
curve represents the postulated experimental path; QCP is 
the corresponding quantum critical point. The Amorphous 
Insulator crosses over from an Anderson-Mott description at 
low Ts to Efros-Shklovskii behaviour at large r^. 

All the transition lines are continuous quantum phase 
transitions for which scaling holds. The point A" is a 
tetracritical point, while is the first order transition of 
the pure two-dimensional electron gas to a Wigner crys- 
tal at which scaling does not hold. In a Coulomb system, 
the density is determined by the charge-neutrality con- 
dition. Hence, the density is continuous; only the order 
parameter is discontinuous. The spin-liquid is a phase 
within the Wigner glass. The Fermi liquid behavior of the 
D = line is unstable with respect to arbitrarily weak 
disorder. It either flows towards the Mott-Anderson in- 
sulating state or to a non-Fermi liquid state. It is thus 
possible that the Mott-Anderson insulator to non-Fermi 
liquid phase boundary terminates in the D ^ Q limit at 
the origin {D = Ts =0), or it may terninate at a point 
of finite interaction strength {D = 0, rg = rf). Also, 
note the reentrance of the Mott-Anderson insulator as a 
function of r^. 

The Wigner glass phase is stable to small thermal per- 
turbations and must exhibit a finite temperature phase 
transition with a transition temperature that vanishes at 
the quantum critical point. However, neither the Mott- 
Anderson insulator, nor the non- Fermi liquid are distinct 
phases of matter at finite temperatures. Because of the 
continuous Heisenberg symmetry of the spin Hamilto- 



nian, there cannot be any magnetic long range order 
in d = 2 at any finite temperature. The transitions 
shown in the figure are true thermodynamic phase tran- 
sitions at T = 0, and are, in principle, detectable as non- 
analyticities in thermodynamic functions, such as the 
chemical potential. A prediction is that the compress- 
ibility, K, of the electrons will vanish as k ^ (n — ricY , 
which should be experimentally observable in capaci- 
tive measurements. (As defined earlier, ly is the cor- 
relation length exponent that has been determined ex- 
perimentally to be 1.5 from the temperature and the 
electric field scalings Q.) The argument is quite sim- 
ple. From macroscopic electrodynamics, the dielectric 
function, e{q,ui = 0) = 1 -|- ^^^K{q). In the insulating 
Wigner glass state, the q = dielectric function is finite 
and different from that of the vacuum. This implies that 
K{q) ~ g as 9 — > 0. Since K{q, ^) should be continuous 
across g^ 1 , it follows from the scaling hypothesis that 
K{q,^) ^ |- in the critical regime, g^ > 1. The critical 
behaviour of the compressibility should be observable in 
macroscopic measurements on the conducting side of the 
transition. 

The magnetic transitions will be reflected in the mag- 
netic response of the electron spins as T ^ 0. In partic- 
ular, a signature of the the spin-liquid to ferromagnetic 
phase transition can be detected in the uniform suscepti- 
bility, which will diverge as the transition is approached 
from the spin-liquid side. A similar signature may be ob- 
servable in the uniform susceptibility at the non-Fermi- 
liquid to Mott-Anderson insulator transition. Because of 
the power-law crystalline order of the Wigner glass, its 
presence can also be detected in surface acoustic wave 
and narrow band noise measurements. 

By applying electric field, a Wigner glass can be made 
to slide |3^. In general, there should not be any asso- 
ciated threshold field; although the equilibrium Wigner 
glass may not contain dislocations, these will be gener- 
ated in the presence of an external electric field. The 
mechanism which destroys any possible threshold field 
is similar to that discussed by Coppersmith js^. The 
current-voltage characteristic in the glass phase is J ~ 
E cxp{—{Ec/ E)'^), where E is the electric field, J is the 
current density, Ec is the typical pinning field, and fi 
is an exponent. This refiects the presence of infinitely 
many low lying metastable states. As the glass melts, 
the current voltage characteristic should become linear 
at asymptotically low currents. Close to the transition, 
the scaling picture suggests that the conductivity can be 
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where 



written as (j{u:,E,T) = a (j^tatt, ^ 
is precisely the exponent x defined earlier. The glass 
correlation length is given by ~ <5~'', and z is the dy- 
namic critical exponent defined by ^ i,g, where i^,- is 
the correlation time. 

We emphasize, once again, that for purer systems, the 
metal-insulator transition must approach the first order 
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melting of a Wigner crystal at which scaling cannot hold, 
while in more disordered samples scaling will hold. The 
scaling theory of the metal insulator transition in- 
volves a transition between an insulator and a perfect 
metal. This description is valid in the critical regime, 
even if the actual asymptotic behavior of the "metallic" 
phase is different from that of a perfect metal; outside the 
critical regime, the behavior of the "metallic" phase can 
be affected by couplings that are irrelevant at the critical 
fixed point (and so do not affect the critical phenomena), 
but are relevant at the putative stable fixed point. Thus, 
it is perfectly consistent with the scaling theory for the 
resistivity outside the critical regime to saturate at low 
temperatures, as has been observed in some experiments, 
or even to turn around and ultimately diverge at very low 
temperatures. 

We hope that many aspects of our proposed phase dia- 
gram can be explored, both theoretically and experimen- 
tally, in the near future. It is also worthwhile to consider 
the extension of the present ideas to three dimensional 
systems, where the existence of a Wigner glass phase is 
on even firmer theoretical footing. 

Note added: We thank V. M. Pudalov for drawing our 
attention to earlier evidence of a Wigner crystal transi- 
tion in Si MOSFETs described in V. M. Pudalov, et at, 
Phys. Rev. Lett. 70, 1866 (1993). However, from the 
perspective of the present paper, this transition is better 
described as a Wigner glass transition with no threshold 
field. 

After this paper was submitted for publication, further 
evidence of a Wigner glass transition was discovered (J. 
Yoon et al, |cond-mat / 9807235| ) . The authors describe 
this as a Wigner crystal transition from the proximity of 
the transition to = 37 obtained from fixed-node quan- 
tum Monte Carlo approximation. This interpretation is 
incorrect as the disorder is finite, which by all theoretical 
accounts must destroy the Wigner crystal. We believe 
that they have observed a transition to the Wigner glass 
state with a power-law order as described here, close, 
however, to the first order melting of the pure crystal. It 
is therefore not surprising that scaling is approximately 
violated; asymptotically close to the transition scaling 
must of course hold as the system is strictly speaking 
not pure. 

H.-W. Jiang in private communications has reported 
to us that the compressibility in his hole doped GaAs 
samples indeed vanishes as a power-law at ~ 8 as 
described above. 

The evidence of the transition to an insulating state 
at smaller values of as shown in Figure 1 has been 
found (A. R. Hamilton et al, |cond-mat/9808108D . The 
character of this second transition is very different from 
that of the first transition. This implies that there must 
be at least two distinct insulating states, as there cannot 
be two different critical transitions to the same insulating 
state from one conducting state. We suggest that this 



lends further support to our phase diagram. 
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